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Abstract 



We show that the production rates of light-flavoured mesons and baryons in inelastic pp 
and pp collisions can be described by a simple approach used to describe data obtained in e + e~ 
annihilation. Based on the idea of string fragmentation, the approach describes the production 
rates of light-flavoured mesons and baryons originating from fragmentation in terms of the 
spin, the binding energy of the particle, and a strangeness suppression factor. Apart from a 
normalization factor and the additional sea quark contribution in inelastic pp and pp collisions, 
pp, pp and e + e~ data at various centre-of-mass energies are described simultaneously. 



1 Introduction 



The soft processes of the fragmentation of quarks and gluons into hadrons cannot be calculated 
with a perturbative approach and instead currently rely on phenomenological models. The 
most successful are the string [|l| and the cluster || fragmentation models implemented in the 
Monte Carlo programs PYTHIA / JETSET §] and HERWIG @J, respectively However, these 
models require either a large number of free parameters in order to reproduce the measured 
hadron production rates (more dramatic in the case of PYTHIA/ JETSET), or do not give 
a satisfactory description of baryon data in e + e~ annihilation, as in the case of HERWIG (a 
review may be found in Ref. || ||). 

In Ref. a simple approach based on the idea of string fragmentation to describe hadron 
production rates in e + e~ annihilation is proposed. We consider that particle production pro- 
ceeds in two stages, namely quark pair production in the colour string field and successive 
recombination. Quark pair production in the colour string field can be considered as a tun- 
neling process. The probability of producing a qq pair is proportional to exp(— nm 2 q / k), where 
m q is the (constituent) quark mass, and k the string constant. We assume that the probability 
of quarks recombining to a hadron with the mass Mh is proportional to exp(— Eund/T), where 
T is the effective temperature in hadronization, and Eu n d = Mh — J2i m qi the hadron binding 
energy, which can be ascribed to the colour- magnetic hyperfine interaction [] . 

The production rates of light-flavoured mesons and baryons from fragmentation can be 
described as 

AT r^l 1 / \N E bind / -t\ 

<N>=C-— (ls) s -e t , (1) 

where 7 S = exp[— •n{m 2 s — m 2 )/^} is the strangeness suppression factor, N s the number of strange 
quarks contained in the hadron, and J the spin of the hadron. C is an overall normalization 
factor, which increases with the increasing centre-of-mass energy (reflecting the rise of multi- 
plicities with increasing energy, which can be predicted by QCD and Cb is the relative 
normalization factor between mesons and baryons (for mesons Cb = 1)- Our approach describes 
quite well the existing e + e~ data on hadron production at various centre-of-mass energies. Fur- 
thermore, it can be applied to heavy flavour production and its predictions there also agree 
well with data. 

Multiparticle production in hadron-hadron interactions with low-momentum transfer, which 
is still one of the least understood processes in high-energy physics, has been studied and 



compared with particle production in e + e~ annihilation by many authors [TCJ, [11], |T2| [13], [L4 
It is generally believed that the mechanism of hadronization is the same in all high-energy 
processes, and that all differences in particle composition can be traced to different initial 
parton configurations Q. It is therefore important to test our approach with data obtained in 
hadron-hadron collisions. 

In this paper, we first show some similarities between the hadron production rates measured 
in inelastic pp collisions and in e + e~ annihilation. Then, taking into account the different initial 



^^One could consider that the production probability of a hadron with a given quark content similar to the 
distribution of the number of atoms at different energy levels of the hyperfine splitting is determined by the 
Boltzmann distribution. 

2 The Lund string model for example, is one explicit realization of this concept. 
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parton configurations and the contribution from sea quarks, we apply our approach to data 
obtained in inelastic pp and pp collisions. 



2 Data on hadron production in pp and pp collisions 



The data on hadron production in inelastic pp and pp collisions used in this analysis are listed 
in Tables 1-4. In some papers, listed in the Tables, pp data are given in terms of production 
cross sections instead of production rates. To calculate the production rate of hadrons per 
inelastic pp collision, we use the value of the total inelastic pp cross sections quoted in the 



papers, or the value from Ref. at the corresponding centre-of-mass energy if inelastic pp 
cross sections are not given in the papers. When counting the hadron production rates, the 
decay products of K^, A, S ± , H ' - and Q~ (and their antiparticles) are not included in the 
papers. This is different from the counting procedure used for e + e~ data. 

The most complete set of data on hadron production in inelastic pp collisions consists of the 
LEBC-EHS results @ at Plab — 400 GeV (centre-of-mass energy ^/s = 27 A GeV), and the 



Fermilab 30-inch bubble chamber data [17 at p LA B — 405 GeV (y/s = 27.6 GeV). This data set 



can be compared with the LEP data fL8fl , the most complete set of data on hadron production 
in e + e~ annihilation. Although the type of interactions and the initial parton configuration 
in pp interactions with low-momentum transfer are different from those in e + e~ annihilation, 
for a direct comparison we can however use hadrons which do not contain the valence quarks 
(u,d) in the proton, such as us, ds and ss mesons, and antibaryons. In most of the cases, such 
hadrons are either produced from fragmentation, or are decay products of other hadrons which 
are produced from fragmentation On the other hand, the contribution of the primary qq pair 
to the production rate of light-flavoured hadrons is small at LEP. It is therefore possible to use 
the following hadrons to compare directly hadron production from fragmentation in inelastic 
pp collisions and in e + e~ annihilation: 



Hadron n pp (27.5GeV)/n ee (91GeV) 

K~ 0.19 ±0.01 

K*" 0.25 ±0.04 

K*° 0.24 ±0.04 

0.18 ±0.02 



Hadron n pp (27.5GeV)/n ee (91GeV) 
p 0.13 ±0.01 

A 0.11 ±0.02 

A 0.20 ±0.13 

£± 0.35 ±0.13 



For all hadrons except the p and A, the ratio n pp /n ee is about the same within errors, indicating 
that hadron production from fragmentation in inelastic pp collisions and in e + e~ annihilation 
differs only by a normalization factor, i.e., the available energy in fragmentation. For the p 
and A, this ratio is lower than that for other hadrons, mainly due to the different counting 
procedures used for decay products of the A, T, ± , E°~ and Q~ (and their antiparticles) as 
mentioned above, and to the additional contribution to the p and A from decays of c and b 
baryons in e + e~ annihilation. If these two corrections are taken into account, the ratio n pp /n ee 
is then equal to 0.21 ± 0.02 for the p Q. From the values of the ratio n pp /n ee listed above, 

3 The contribution of sea quarks and antiquarks to hadron production is small (see discussions in the next 
section) 

4 For the A the uncertainty of the corrections is large since the branching fractions c(b) baryons — > A are not 
well known. 
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we see that the number of hadrons produced from fragmentation in inelastic pp collisions at 
s/s ~ 27.5 GeV is about 20% of that in e + e~ annihilation at y/s = 91 GeV. We will discuss 
this in more detail in the next section. 



3 Analysis 

We first consider hadron production in inelastic pp collisions. Similar to hadron production 
in e + e~ annihilation, this process can also be considered to proceed in four steps: interaction 
between two initiator partons (valence quarks, gluons, sea quarks or antiquarks) with two beam 
remnants left behind, followed by a parton shower development from the initiator partons, and 
subsequently the transition from partons to hadrons. Finally, the unstable hadrons decay 
according to their branching ratios. 

For a proton the possible scenarios for the initiator parton and beam remnant are the 
following: 

1. If the initiator parton is a valence quark q v , the beam remnant is simply a diquark 
composed of the two leftover valence quarks, i.e. either a ud or a uu diquark. The 
initiator parton and the diquark are at the two end-points of a string which is then 
hadronized in the usual way ||. 

2. If the initiator parton is a sea quark q s or antiquark q s , the beam remnant contains four 
quarks: uudg s or uudg s . Since the q s q s pair, to a first approximation, is in a colour-octet, 
the subdivision uud+g s (or uud+g s ) is not allowed, since it would correspond to a colour- 
singlet q s q s ||. Therefore uud will be subdivided into a quark (u or d) and a diquark 
(ud or uu). In this case in the beam remnant we have a diquark and two single quarks 
(or a quark and a antiquark), which will then be at the end-points of two strings f\ and 
hadronized in the usual way. We introduce two free parameters for the fractions of the 
sea quarks: f = f u = fd and x s = f s /f. 

3. If the initiator parton is a gluon g, the beam remnant is a colour-octet uud state, which 
is subdivided into a quark and a diquark. The treatment is similar to that for scenario 1. 

On average, the individual numbers of the primary quarks and diquarks in an inelastic pp 
collision are: N u = 2(§ + /), N d = 2(| + /), jV s = N- d = 2/, N B = N- B = 2x.f, N uu = 2 x § and 
N uA = 2 x §. 

Taking into account the numbers of the primary quarks and diquarks, hadron production 
in inelastic pp collisions can also be described by our approach discussed in Ref. 0. For 
simplicity, as in Ref. J7J, we apply our approach analytically to the data. Instead of using 
the parton distributions p(x, Q 2 ) (p = q v ,q s ,q s , g), for the above three scenarios we use an 
averaged parameter C in Eq. ([!]), which is integrated over x. The total hadron rates are 
calculated as follows. Firstly we calculate the number of primary hadrons. Primary hadrons 
are defined as those which are not decay products of other hadrons, i.e. as hadrons originating 

5 Here we use the treatment of sea quarks and antiquarks as discussed in Ref. |fl9|| , which is different from 
the conventional way [[| . 
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from fragmentation, or containing a primary quark or a diquark from the p in an inelastic pp 
collision. We use Eq. (HD to calculate the number of light-flavoured hadrons produced from 
fragmentation. Heavy quark production from fragmentation is strongly suppressed due to the 
term exp(— irrri^/ k) and can therefore be neglected. For hadrons which contain a primary quark 
q(q = q v , q s , q s ) or a primary diquark q\q\, we use Eq. (|TJ) to determine their relative ratios, 
and then obtain their rates by normalizing the sum of the rates to the average number of q or 
q\q 2 in an inelastic pp collision. Compared to the production cross section of light-flavoured 
hadrons, the charm and bottom production cross sections in pp collisions are small ]2D|] and 



can therefore be neglected in this analysis. The spin of the diquark is taken into account in 
the calculation. For example, the A contains a spin-0 ud and the S° a spin-1 ud diquark. All 
light-flavoured hadrons up to a mass of 2.5 GeV in the meson and baryon summary table of 
Ref. jn| are included in the calculation. Finally, we let all the primary hadrons decay according 



to their decay channels and branching ratios given in Ref. |T3[ . The decay chain stops when \i 



7i, K , K° L , K° s , A, E*, S 0, , n (or their antiparticles) or stable particles are reached. 

In the fit to pp data at y/s = 27.4-27.6 GeV, we choose C, f and x s as free parameters. For 
the other parameters in Eq. ([!]) we use the result of the fit to e + e~ data [0 : 7s = 0.29 ± 0.02, 
Am = 0.161 ± 0.024 (GeV), T = 0.298 ± 0.015 (GeV) and C B = 11.0 ± 0.9. In the fit 
these parameters are set to their central values since otherwise there would be too many free 
parameters which are highly correlated with each other, leading to unstable fit results. The 
value of m s is set to 0.5 GeV 0. The fit is performed by minimizing the error function 

X 2 = Y,{N- alc ~ A/™) 2 /(AA™) 2 , (2) 

i 

where Nf alc is the calculated rate, j\T meas and AAf meas the measured rate and its error, respec- 
tively. 

The fit results are listed in Table 5 and shown in Fig. 1. The value of the normalization 
factor C is equal to 0.042 ± 0.007, where the error also includes the uncertainty on 7 S , Am, T 
and Cb which are set to their central values in the fit. The ratio of this value and the value 
of C obtained from the fit to e + e~ data at -^i = 91 GeV [7J] is (19 ± 3)% (common systematic 
errors have been taken into account), which is in a good agreement with the value of n pp /n ee 
discussed in the previous section. From the fit, the probability that the initiator parton is a u 
(or d) sea quark is determined to be / = 0.044 ± 0.010. The strange sea content is determined 



to be x s = f s /f = 0.51 ± 0.22, which is consistent with the parametrization used in Ref. |2T 



If we also choose the parameter T as a free parameter in the fit, we obtain T = 0.271 ± 
0.010 GeV, which is in agreement with the values of T obtained from the fits to e + e~ data 0. 
This indicates that the value of T is universal in different types of interactions. 

As can be seen from Table 2 and Fig. 1, for most of the hadrons the calculated rate agrees 
well with the measurement (within two standard deviations). Only at three data points, namely 
7r + , p° and is the difference between the calculated and measured rates larger than three 
standard deviations. These three data points contribute more than half of the \ 2 value of 
the fit (49 out of the total \ 2 °f 85 for 33 data points, fitted with three free parameters). 
Another major reason for the large \ 2 value is because the uncertainty of the calculated rates 

6 Since exp(-E bind /T) = exp(J2i m u /T) ■ exp{-[M h - J2ii m n ~ m u)]/T}, the factor exp(^ 4 m u /T) can be 
absorbed in C and Cb- The error function of the fit is mainly sensitive to the change in the mass difference 
m „ - to,,,. 
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due to uncertainties on the mass, branching ratios and decay modes of the resonance states is 
not taken into account in the fit. If, for example, we introduce a 5% error on the number of 
secondary particles as the uncertainty of the calculated rates (AN£ alc = 0.05* (1 — /f" m ) *N^ alc , 
where f? rim is the primary fraction of hadron i) and add it to the term AiV™ s in Eq. (0), 
the x 2 /d°f value of the fit is then reduced from 85/30 to 61/30, while the fit results remain 
essentially unchanged. Comparing the fit to almost the same data set with the same number 
of free parameters in Ref. fL4"l , the \ 2 /dof value of our fit is much lower [] and our calculated 



rates are closer to the measured rates. 

In the fit to pp and pp data at the other centre-of-mass energies we choose only C as a 
free parameter because of the small number of data points available at some energies. The 
parameters / and x s are set at the values obtained from the fit to pp data at y/s = 27.4- 
27.6 GeV. In the pp data set at y/s = 52.5-53.0 GeV, f 2 (1270), K^(1430) and Ka ± (1430) are 



not included in the fit as the results from Ref. |22[ and |23| differ significantly. The treatment 



of initiator partons and beam remnants in pp collisions is similar to that in pp collisions, with 
one of the two initiator partons and beam remnants replaced by their antiparticles. The fit 
results are listed in Tables 1, 3, 4 and 5. One can see that the calculated rates agree well with 
the measurements. 

The fraction of primary hadrons is also given in Tables 1-4. Owning to different initial 
parton configurations, the fraction values in Tables 1-4 are different from the corresponding 
values for hadrons produced in e + e~ annihilation at y/s = 91 GeV (see the last column in 
Table 2): 



• For mesons, and in inelastic pp collisions also antibaryons, the fraction value for pp (pp) 
data is in general slightly higher than that for e + e~ data, mainly due to the fact that 
the amount of c and b hadrons produced in pp (pp) collisions is negligible. In contrast, 
in e + e~ annihilation at yfs = 91 GeV about 40% of the events contain c or b hadrons 
which decay to light-flavoured hadrons, leading to a lower value of the primary fraction 
of light-flavoured hadrons. 

• For baryons which contain a uu (uu) or ud (ud) diquark, the fraction value for pp (pp) 
data is much higher than that for e + e~ data as a uu (uu) or ud (ud) diquark already 
exists in the initial configuration of inelastic pp (pp) collisions. In contrast, for ddd (ddd) 
or dds (dds) types of baryons, for example, the £~, the fraction value can be lower than 
that for e + e~ data as such baryons cannot be produced directly from the diquark in the 
p (p), but can be decay products of the other baryons which are produced directly from 
the diquark in the p (p). 



4 Conclusions 

We have shown that the production rates of light-flavoured mesons and baryons in inelastic pp 
and pp collisions are described by a simple approach used to describe e + e~ data. Taking into 
account the different initial parton configurations and the additional sea quark contribution in 

7 Note in Ref. (l4) a AN£ alc term has been taken into account in the fit. fn some cases, for example, for p° 
and A++, ANf alc is much larger than AN™ eas . 
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inelastic pp and pp collisions, pp, pp and e + e~ data on the production of light-flavoured hadrons 
at various centre-of-mass energies are described simultaneously (apart from a normalization 
factor, which reflects the rise of multiplicities with increasing energy). Moreover, as shown 
in Ref. H, data on heavy flavour production in e + e~ annihilation are also described by our 
approach. All this shows that our approach can provide a universal description of hadron 
production, irrespective of the type of the interaction and the initial parton configuration in 
various interactions. 
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Hadron 



Rate 
Measured 



Rate 
Calculated 



Primary Fraction 
Calculated 



References 



y/s = 19.4-19.7 GeV {o mel = 32.1 mb p4 



Charged 
Neg. charged 

7 

P° 
A 

A 



7.68 ±0.06 7.53 

2.85 ±0.03 2.90 

6.68 ±0.48 6.75 

0.174 ±0.011 0.205 

0.33 ±0.06 0.39 

0.098 ±0.010 0.124 

0.014 ±0.004 0.011 



0.28/0.30^ 
0.49 
0.27 
0.15 



ES El 
11 H 



2| 



y/s = 23.3-23.7 GeV (a mel = 32.21 mb |28 



Charj 

7T° 

7T + 

K+ 

K" 
K *± 

P 
P 
A 

A 

s *± 



;ed 9.24 ± 1.39 8.61 

3.42 ±0.62 3.79 0.17 

3.71 ±0.37 3.67 0.19 

3.27 ±0.33 3.14 0.17 

0.214 ±0.025 0.250 0.28/0.30^ 

0.337 ±0.051 0.344 0.33 

0.209 ±0.031 0.212 0.28 

0.137 ±0.043 0.215 0.60/0.67^ 

1.63 ±0.24 1.10 0.24 

0.085 ±0.013 0.063 0.16 

0.112 ±0.016 0.131 0.27 

0.020 ±0.004 0.018 0.15 

0.017 ±0.012 0.021 1.00 

0.014 ±0.011 0.004 1.00 



1 

in 



mi 

m 

in 







y/s = 26.0 GeV (a inel = 32.8 mb |32j) 



Charged 
Neg. charged 

A 
A 



9.06 ±0.09 8.92 

3.53 ±0.05 3.60 

0.26 ±0.01 0.26 

0.12 ±0.02 0.13 

0.013 ±0.004 0.020 



0.28/0.30^ 
0.26 
0.15 



(a) Primary fraction of K° and K U , respectively. 

(b) Primary fraction of K*~ and K* + , respectively. 



3 
PI 



m 

in 

m 



Table 1: Average hadron production rates per inelastic pp collision at various centre-of-mass energies 
(excluding charge conjugates and antiparticles if not indicated), compared with the calculated values. 
The fraction of primary hadrons obtained from the ht is also shown. For the calculated rates and 
fractions see discussions in Section 3. 
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Hadron Rate 


Rate 


Prim. Fraction References 


Prim. Fraction 


Measured 


Calculated 


Calculated 


e+e" 


v/i = 27.4-27.6 GeV (a mel 


= 32.8 mb m; 







7T 
7T + 

7T~ 

K+ 

K 

r/ 

P° 

P + 

P 

K*° 
K*° 
K*+ 

K*- 



fo(980) 

f 2 (1270) 

K^(1430) 

P 

P 
A 

A 

£+ 

s- 

A++ 
A— 
A 
A 

£*+ 

g*± 

A(1520) 



3.87 ±0.12 3.82 

4.10 ±0.11 3.70 

3.34 ±0.08 3.17 

0.232 ±0.011 0.252 

0.331 ±0.016 0.346 

0.224 ±0.011 0.215 

0.30 ±0.02 0.33 

0.385 ±0.018 0.460 

0.552 ±0.082 0.490 

0.354 ±0.058 0.408 

0.120 ±0.021 0.117 

0.090 ±0.016 0.079 

0.132 ±0.016 0.137 

0.088 ±0.012 0.080 

0.391 ±0.024 0.421 

0.0189 ±0.0018 0.0175 

0.023 ±0.008 0.035 

0.092 ±0.012 0.090 

0.11 ±0.05 0.04 

1.20 ±0.097 1.11 

0.063 ±0.002 0.064 

0.125 ±0.008 0.131 

0.020 ±0.004 0.018 

0.048 ±0.015 0.039 

0.013 ±0.006 0.021 

0.218 ±0.0031 0.214 

0.0128 ±0.0049 0.0105 

0.141 ±0.008 0.134 

0.0336 ±0.010 0.011 

0.0204 ±0.0024 0.0193 

0.0100 ±0.0018 0.0022 

0.0078 ± 0.0025 0.0045 

0.0171 ±0.0031 0.0161 



0.28 



0.17 
0.19 
0.17 
./0.30< a > 
0.33 
0.28 
0.28 
0.50 
0.54 
0.48 
0.59 
0.59 
0.67 
0.60 
0.52 
1.00 
1.00 
0.77 
1.00 
0.24 
0.16 
0.27 
0.15 
0.48 
0.10 
0.83 
0.72 
0.73 
0.71 
1.00 
1.00 
1.00 
0.83 



m 











in 
in 
in 
in 

in 
in 

16 



16 

in 

17 



16 





in 
in 


m 



16 







TT 



= 

n 

16 


in 




0.16 
0.18 
0.18 
0.25 
0.23 
0.23 
0.30 
0.47 
0.58 
0.58 
0.49 
0.49 
0.50 
0.50 
0.48 
0.64 
0.99 
0.78 
1.00 
0.12 
0.12 
0.12 
0.12 
0.42 
0.42 
0.69 
0.69 
0.69 
0.69 
0.91 
0.91 
0.91 
0.71 



(a) Primary fraction of K and K , respectively. 



Table 2: Average hadron production rates per inelastic pp collision at yfs = 27.4—27.6 GeV (excluding 
charge conjugates and antiparticles if not indicated), compared with the calculated values. The fraction 
of primary hadrons obtained from the fit is also shown. For the calculated rates and fractions see 
discussions in Section 3. For comparison, the fraction of primary hadrons obtained from the fit to 
e + e~ data at 91 GeV is also given in the last column (the fraction values have been averaged for 
baryons belonging to the same isospin multiplet). 
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Hadron 



Rate 
Measured 



Rate 
Calculated 



Primary Fraction 
Calculated 



References 



y/s = 30.6 GeV 


Charged 


10.07 ± 1.51 


9.87 


7T + 


4.07 ±0.41 


4.21 


7T~ 


3.65 ±0.37 


3.68 


K+ 


0.367 ±0.055 


0.400 


K" 


0.244 ±0.037 


0.269 


P 


1.63 ±0.24 


1.13 


P 


0.108 ±0.016 


0.091 



0.20 
0.18 
0.33 
0.29 
0.24 
0.16 



m 

HI 

m 



m 

HI 

hi 

HI 



y/s = 44.6 GeV 



Charged 10.99 ± 1.65 10.82 

7T+ 4.45 ± 0.45 4.62 

7T 4.09 ±0.41 4.09 

K+ 0.411 ±0.062 0.443 

K" 0.286 ±0.043 0.312 

p 1.62 ±0.24 1.15 

p 0.132 ±0.020 0.112 



0.20 
0.18 
0.33 
0.29 
0.24 
0.16 



m 
m 



m 

Hi 
Ell 
HI 

91 



y/s = 52.5-53.0 GeV (a inel = 35.0 mb |T3 



Charged 


11.47 ± 1.72 


11.24 




7T + 


4.68 ±0.47 


4.80 


0.20 




4.29 ±0.43 


4.27 


0.18 




0.329 ±0.071 


0.359 


0.29/0.31 


K+ 


0.430 ±0.065 


0.462 


0.33 


K~ 


0.306 ±0.046 


0.330 


0.29 


P° 


0.63 ±0.14 


0.62 


0.51 


K*° 


0.123 ±0.028 


0.123 


0.60 




0.037 ±0.010 


0.026 


1.00 


f 2 (1270)W 


0.154 ±0.034 


0.122 


0.77 


f 2 (1270)W 


0.075 ±0.007 


0.122 


0.77 


K^(1430)^ 


0.0044 ±0.0017 


0.0264 


1.00 


Ka ± (1430) W 


0.028 ±0.007 


0.0205 


1.00 


Ka ± (1430) W 


0.0031 ±0.0015 


0.0205 


1.00 


P 


1.62 ±0.24 


1.16 


0.24 


P 


0.144 ±0.022 


0.121 


0.16 



HI 
HI 



HI 

m\ 

HI 
HI 
HI 
HI 
HI 
HI 



HI 



HI 
HI 
HI 



HI 



(a) Primary fraction of K° and K 

(b) Not included in the fit. 



respectively. 



Table 3: Average hadron production rates per inelastic pp collision at various centre-of-mass energies 
(excluding charge conjugates and antiparticles if not indicated), compared with the calculated values. 
The fraction of primary hadrons obtained from the ht is also shown. For the calculated rates and 
fractions see discussions in Section 3. 
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Hadron Rate 


Rate 


Primary Fraction References 


Measured 


Calculated 


Calculated 


y/s = 200 GeV 



H 1 ^ 



Charged 



n 
A 



21.4 ±0.4 21.4 

0.75 ±0.09 0.78 

0.75 ±0.10 0.73 

0.23 ±0.06 0.16 

0.015 ±0.015 0.007 



0.31 
0.19 
0.20 
0.50 



3| W 



m 

13 



= 546 GeV 



Charged 
A 



29.4 ±0.3 29.4 

1.12 ±0.08 1.12 

0.265 ±0.055 0.205 

0.05 ±0.015 0.011 



0.31 
0.19 
0.50 



34 w 



3j 



H 



y/S = 900 GeV 



Charged 



n 
A 



35.6 ±0.9 35.7 

1.37±0.13 1.38 

1.0 ±0.2 1.1 

0.38 ±0.08 0.24 

0.035 ±0.020 0.014 



0.31 
0.18 
0.18 
0.50 



3J 
3J(6) 



3j 



3J] 



(a) The average charged multiplicity value quoted in this reference 

is increased by one to include the leading particles. 

(b) The average production rate of the neutron quoted in this reference 

is increased by 0.5 to include the leading particles. 



Table 4: Average hadron production rates per non-single-diffractive pp event at various centre-of- 
mass energies (excluding charge conjugates and antiparticles) , compared with the calculated values. 
The fraction of primary hadrons obtained from the fit is also shown. For the calculated rates and 
fractions see discussions in Section 3. 
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v^(GeV) 


C 


/ 


X s 


X 2 /dof 


PP 


19.4-19.7 


0.026 ± 0.004 


0.044 (fixed) 


0.51 (fixed) 


25/6 


23.3-23.7 


0.041 ±0.006 


0.044 (fixed) 


0.51 (fixed) 


16/13 


26.0 


0.046 ±0.007 


0.044 (fixed) 


0.51 (fixed) 


7/4 


27.4-27.6 


0.042 ±0.007 


0.044 ±0.010 


0.51 ±0.22 


85/30 


30.6 


0.060 ±0.010 


0.044 (fixed) 


0.51 (fixed) 


6/6 


44.6 


0.074 ±0.014 


0.044 (fixed) 


0.51 (fixed) 


6/6 


52.5-53.0 


0.080 ±0.010 


0.044 (fixed) 


0.51 (fixed) 


7/10 


PP 


200 


0.229 ±0.035 


0.044 (fixed) 


0.51 (fixed) 


2/4 


546 


0.347 ±0.053 


0.044 (fixed) 


0.51 (fixed) 


8/3 


900 


0.439 ± 0.066 


0.044 (fixed) 


0.51 (fixed) 


4/4 



Table 5: Results of the fit to pp and pp data obtained at various centre-of-mass energies. The errors 
given in the Table also include the uncertainty on the parameters which are set to their central value 
in the fit. 
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Figure 1: The upper plot shows the average hadron production rates of light-flavoured hadrons 
originating from fragmentation at y/s = 27.4-27.6 GeV in inelastic pp collisions (measured value x 
fraction of hadrons originating from fragmentation as determined by the fit), multiplied by the factor 
Cb/[7^ Vs (2J + 1)] (see Eq. (|l|)), as a function of the binding energy of hadrons. The fit results are 
shown as the line. The lower plot shows the difference between the predicted and measured rates in 
terms of the number of standard deviations. 
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